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tin, ABSTRACT 

' Direct and indirect observational evidence leads to the conclusion that high redshift 

QSOs did shine in the core of early type proto-galaxies during their main episode 
' of star formation. Exploting this fact, we derive the rate of formation of this kind of 

\ stellar systems at high redshift by using the QSO Luminosity Function. The elemental 

'b) . proportions in elliptical galaxies, the descendents of the QSO hosts, suggest that the 

^—^ ' star formation was more rapid in more massive objects. We show that this is expected 

" , \ to occur in Dark Matter haloes, when the processes of cooling and heating is con- 

sidered. This is also confirmed by comparing the observed sub-mm counts to those 
' derived by coupling the formation rate and the star formation rate of the spheroidal 

Q>^ , galaxies with a detailed model for their SED evolution. In this scenario SCUBA galax- 

ies and Lyman Break Galaxies are early type proto-galaxies forming the bulk of their 
stars before the onset of QSO activity. 
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• '-j ' 1 INTRODUCTION sion and absorption lines show that high redshift QSOs live 
Itn] in a metal enriched environment (see Hamann & Ferland 
H ; QSOs had been for a long time the main probe of the ^ggg ^ comprehensive review). Observations at far-IR 
. 5y , epoch when galaxies are thought to have formed. Hubble ^nd sub-mm wavelengths uncovered large amounts of dust 
Space Telescope observations, especially the HDF surveys, ggOs (Omont et al. 1996; Ivison et al. 1998; Benford et al. 
and ground based observations with the new 10 meter class ^ggg. yun et al. 2000). The dust emission has been ascribed 
telescopes, have opened to direct investigation also the early ^o starburst in the host galaxy (see e.g. Yun et al. 2000; 
phases ofgalaxy formation. Exploration of the local Universe Chapman et al. 1998; Benford et al. 1999) or dust iUumi- 
is also yielding extremely relevant results to understand ^^^^^ by the active nucleus (see e.g. Andreani, Franceschini 
the QSO phenomenon. Very high angular resolution pho- ^ Granato 1999). The evidence of the dust illuminated by 
tometric and spectroscopic observations have demonstrated starbursts in high redshift radio galaxies is solid (see e.g. 
that Massive Dark Objects (MDO) are generally present in Archibald et al 2000) 
nearby galaxies endowed with significant spheroidal compo- 
nents (Magorrian et al. 1998; van der Marel 1999). MDOs ^.11 these pieces of evidence suggest that QSOs at high 
are thought to be dormant BHs which spent their shin- redshift are active nuclei shining in early type galaxies, dur- 
ing phase as QSOs. Indeed, the MDO mass function well i^ig ^he short fraction of the Hubble time when they were 
matches the mass function of baryons accreted onto BHs vigorously forming stars and still gas rich. The high metal 
during the QSO activity (Salucci et al. 1999, hereafter paper abundance of the QSO environments strongly supports the 
I) . Estimated MDO masses are roughly proportional to those i<iea that the bulk of star formation in host galaxies occurred 
of the spheroidal component of the host galaxy (Magorrian before the QSO shmmg phase. 

et al. 1998; Gebhardt et al. 2000) . This fact imphes that in fact there is evidence that the bulk of stars have 

high redshift and highly luminous QSOs have been hosted formed very soon in elliptical galaxies and in the prominent 

in massive early type galaxies. bulges of the spirals. Cluster elliptical galaxies exhibit a very 

In addition, spectroscopic observations of QSO emis- tight color-magnitude relation, which can be explained if the 
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bulk of their stars formed at very early epochs, correspond- 
ing to z >2 (Bower et al. 1992; Ellis et al. 1997; Kodama et 

al. 1998). The very weak dependence of the Mg2 — (J relation 
on the galaxy environment suggests that field ellipticals are 
on average at most 1 Gyr younger than cluster ellipticals 
(Bernardi et al. 1998). The bulges of spiral galaxies exhibit 
a tight correlation between the Mg2 index and the intrinsic 
luminosity similar to that of the ellipticals (Jablonka, Mar- 
tin & Arimoto 1996). These facts support the conclusion 
that spheroids formed stars very rapidly at early epochs (see 
e.g. Renzini 1999). The rapid star formation can be strongly 
reduced or halted by the exhaustion of the gas in most mas- 
sive objects or by heating and winds, if SN explosions (or the 
QSO activity itself, as argued below) transfer to ISM even 
a fraction of their energy (see e.g. Dekel and Silk 1986). For 
instance galactic winds were brought in the play in order 
to explain the mass-metallicity relation in elliptical galaxies 
( Mathews & Baker 1971; Larson 1974). The observed iron 
abundances in groups and clusters of galaxies can be ex- 
plained if galactic winds expelled large amounts of enriched 
gas from spheroids (see e.g. Renzini 1997). Many photomet- 
ric studies confirm that the massive early type galaxies were 
practically formed at z ~ 1 with a number density very 
close to the local one (Im et al. 1996; Franceschini et al. 
1998; Schade et al. 1999), even though the issue is still some- 
what controversial (e.g. Kauffmann, Chariot & White, 1996; 
Fontana et al. 1999). 

In this paper we present evidence in favour of a strict re- 
lationship in the formation of QSOs and their spheroid hosts, 
with particular emphasis on the timing of star formation in 
the hosts and the QSO shining. We derive the formation rate 
of the spheroids from the formation rate of the QSOs (sect 
2.1 and 2.2), taking into account the observational and the- 
oretical supports for higher star formation rates (SFRs) in 
more massive spheroids. In section 3 we present the model 
(GRASIL) we use to estimate the evolution of the spectral 
energy distribution (SED) of the spheroids. The model in- 
cludes dust. In section 4, exploting the formation rate of 
spheroids coupled with GRASIL, we argue that the QSO 
hosts at z > 2 are the galaxies uncovered in the submillime- 
tre surveys with SCUBA (Blain et al. 1999a, b,c). In section 
5 we discuss relationship of the young spheroids with Lyman 
Break Galaxies (Steidel, Pettini and Hamilton 1995; Steidel 
et al. 1996; Madau et al. 1996) and in section 6 we investi- 
gate the QSO phase. In section 7 we place our results in the 
context of hierarchical cosmogonies. Section 8 is devoted to 
a summary. 

Unless otherwise specified, the results we present have 
been obtained adopting a cosmological model with Hq = 
70 km s"^ Mpc"\ = 0.7, and Q.m = 0.3. We also per- 
formed the same computations in an Einstein-de Sitter cos- 
mology with Ho — 50. The results, being essentially iden- 
tical, are not shown. To avoid confusion, we indicate with 
lowercase t, times measured from the Big Bang, while galac- 
tic ages (i.e. times measured from the onset of star formation 
in a galaxy) are indicated with uppercase T. 



2 THE FORMATION OF SPHEROIDAL 
GALAXIES AS TRACED BY QSO 
EVOLUTION 

The epoch-dependent QSO Luminosity Function (LF) 
nQso{LQso, z) is well defined in a large redshift interval 
< z < 3 both in the optical and in the soft X-ray bands 
(Boyle et al. 1993; Pei 1995; Kennefick et al. 1996; Grazian 
et al. 2000; Miyaji, Hasinger and Schmidt 2000); also new 
optical surveys allow sound estimates of the QSO evolu- 
tion at higher redshift. (Fan et al. 2000). Observations with 
Chandra X-ray Observatory have shown that most of the 
Hard X-ray Background (HXRB) is due to active nuclei 
( Mushotzky et al 1999; Giacconi et al 2000; Barger et al 
2000). In the following section we will use the Luminosity 
Function of QSOs, including the 'obscured ' ones responsi- 
ble for a large fraction of the HXRB, in order to infer the 
formation rate of the host spheroidal galaxies. 

2.1 The formation rate of spheroids 

Since quasar lifetimes Atq are short compared to the typical 
evolution timescale of their LFs, the rate at which QSOs 
with BH mass Al, shine at time tgso (corresponding to 
redshift zqso) is related to the LF by: 



h.{M.{LQso),tQso) 



nQsoiLqso, zqso) diqso 



AtQ 



dM. 



(1) 



Here nqso is taken from Pei (1995) for the optical bands, 
while the contribution from obscured objects is estimated 
using the LF of Boyle et al. (1993) (see Comastri et al. 1995 
and paper I for details). The luminosity of a QSO in a given 
e.m. band is related to the central BH mass Al, through: 



LqS0{M.) = /EDLED(Af.)/Cs, 



(2) 



where Cs is the bolometric correction appropriate for the 
e.m. band used and /ed = iboi/^ED is the bolometric lumi- 
nosity in units of the Eddington luminosity. The available 
data on AGN and QSOs suggest that the ratio /ed is a 
function of the redshift and/or of the luminosity, going from 
- 0.05 - 0.1 for faint local AGNs, to > 1 for bright high 

redshift QSOs (see, e.g., Padovani 1989; Sun and Malkan 
1989; Wandel 1999). This behaviour can be represented as 
function of the luminosity by: 



/ed 



Lhoi 



1049erg/s 



(3) 



with the exponent qed set to 0.2 (cfr paper I). 
We further assume that quasar lifetimes are proportional to 
tduty, the e-folding time for the exponential growth of BH 
mass during QSO activity: 



AtQ = Ktduty = Ke 



(4) 



I/ed/ed 

If we adopt e = 0.1 for the mass to radiation conversion effi- 
ciency of accretion, we have fduty = 4 x 10^ yr for /_bd ~ 1. 
The value Atg can be inferred by matching the local BH 
Mass Function with the Mass Function inferred from ac- 
cretion (Salucci et al 1999), by the QSO clustering proper- 
ties (Martini and Weinberg 2000) or by matching the space 
density of high redshift QSOs and the space density of lo- 
cal massive spheroids (Richstone et al 1998; Monaco et al 
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Figure 1. Shining rate of bulges dn^Msjyh/d-^ redshift 1 (dot- 
ted tliin line), 2 (long-dashed), 3 (dot-dashed) and 4 (short-long 
dashed). The thick continuous line shows the total mass function 
of bulges, n{Mgph), compared to the one inferred in paper I from 
the local LF of galaxies (thick dashed line). 



over time yields their present day mass function. To com- 
pare it with data, we use the following relationship between 
the visible mass of a spheroid Msph and its light: 

Msph/i = {M/L),{L/L,Ef''^''-\ (7) 

where L,e is the Schechter parameter for the E luminosity 
function, P^pb = 1-25 and {M/L)t = 4.9 for luminosities in 
the B band (cfr. paper I). With this recipe the local LF of 
spheroids is translated into a Mass Function (see e.g. Salucci 
et aL 1999). 

We are mainly interested to high redshift, since we want 
to test the hypothesis that the QSO activity signals the end 
of the main episode of star formation in spheroids; in our 
view the low redshift z < 1 activity is mostly due to reac- 
tivation (see e.g. Kauffman and Haehnelt 2000; Cavaliere 
and Vittorini 2000). The galaxies associated with QSOs 
are bright objects. Following the correlation between the 
bulge luminosity and the BH mass (Magorrian et al 1998) 
or the correlation between velocity dispersion and BH mass 
(Ferrarese and Merritt 2000; Gebhardt et al 2000 ), QSOs, 
which typically have Ltot 3 x 10*^ and can be associated 

to massive BH {Mbh ^ 2 x 10^), are expected to have been 
hosted in spheroidal galaxies with central velocity dispersion 
a ^ 100 km/s and with present luminosity Mb ~ — 18. 



1999). For instance the analyses of Martini and Weinberg 
(2000), Salucci et al (1999) and Monaco et al (2000) suggest 
AtQ ~ 8 X 10^ - 10** yr. 

Using equations (^) and (^) into equation (^) we get: 

h,{M,,tQso) = 

nQso(LQso(M.),.gso) ^'^^°^'^-\^ ' , -^(5) 

M, (1 — qed ) Atq 

In order to proceed we introduce here a major assumption, 
namely that the relationship between the BH mass and the 
host mass (or velocity dispersion), holding for local spheroids 
(Magorrian et al 1998; Ferrarese and Merritt 2000; Gebhardt 
et al 2000), has been imprinted during the early phase of 
the QSO and host evolution (see e.g. Silk and Rees 1998; 
Fabian 1999), and has not been changed significantly by 
subsequent mergers. Then h, {M,,tQso) is linked to the rate 
nsph{Msph,tQso) at which galaxies endowed with spheroid 
of present day mass Msph appear to host a QSO at time 
tgso by, 

h.{M.,tQso)^ J nsph{Msph,tQso)f{x)dx, (6) 

where f{x) is the distribution of a; = log{ Data on 

local objects suggest that f[x) can be represented by a 
gaussian with mean (x) = —2.6 and a = 0.3 (see paper I). 
Equating the second members of the two equations above, 
we get a relationship which, once deconvolved, yield an esti- 
mate of hsph{Msph,tQSo), except for the somewhat uncer- 
tain factor l/{Ke). This residual uncertainty can be elim- 
inated normalizing the result in order to recover the local 
LF of spheroids. Fig. |l| shows that there is a good agree- 
ment between predictions and observations, provided that 
Ke ~ 0.1. In the figure the rate hsph{Msph,tQso) at dif- 
ferent redshift is presented. Integration of nsph{Msph,tQso) 



2.2 Early star formation in spheroidal galaxies 

In the Introduction we have reviewed several pieces of evi- 
dence suggesting that star formation begins in the hosting 
spheroids at a time and proceeds vigoursly at least un- 
til the time tgso, when the QSO shines. Here we estimate 
the duration of the star formation phase and its possible 
dependence on the mass. 

The statistics on Broad Emission Lines suggest that 
metallicities are possibly higher in more luminous QSOs 
(Hamann and Ferland 1993, 1999). This is reminiscent of 
the mass-metallicity relationship in the local ellipticals and 
can be naturally explained taking into account the correla- 
tion between QSOs and host galaxy luminosities expected 
on the basis the observed correlation between BH remnants 
and spheroidal masses. 

Large elliptical galaxies exhibit, at least in the central 
regions, overabundances of magnesium and possibly of other 
Q-elements. Fisher et al (1995) have demonstrated that ellip- 
ticals with larger central velocity dispersion ctq tend to have 
larger Mg/Fe. More recently Trager et al (2000b) showed 
that there is a clear correlation between the ratio E/Fe Q 
within r — re/8 and the velocity dispersion ct in a significant 
sample of local ellipticals. Their analysis, extended also to 
other properties, favours a scenario in which local ellipticals 
are composed by an old population encompassing ~ 80% of 
the stars and by a 20% of 'young' stars. In order to explain 
the E/Fe-CT correlation and the Z-plane (linking [Z/H], a 
and age), two possibilities have been suggested: either the 



* "E" refers to the total mass fraction of the elements whose 
abundance is enhanced in giant elliptical galaxies with respect 
to Fe; this group is similar but not exactly the same as the a- 
elements group (for more details see Trager et al 2000a). 
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IMF flattens with increasing a, or a elements are more ef- 
fectively ejected in smaller galaxies due to stronger early 
winds. The two possibilities are not mutually excluding. 

By converse, several authors suggested that in order to 
reproduce the [E/Fe]-cr correlation with reasonable parame- 
ters for the IMF and for the SNIa rate and yields, the time 
scale for the burst must be short (~0.5-l Gyr) for larger 
masses (see e.g. Matteucci 1994; Thomas, Greggio and Ben- 
der 1999). 

In the following we will test a scenario in which the 
star formation begins at and proceeds very rapidly in 
most massive ellipticals and, at tqso, the QSO stops the 
star formation by heating and ionizinig the ISM and/or by 
triggering winds. Afterwards, only minor episodes of star 
formation can occur, in the sense that they involve ^ 20% 
of the mass. In order to obtain the observed dependence of 
E/Fe and Z/H on mass the star formation time scales are 
shorter (or conversely SFRs are larger) for more massive 
galaxies. 

We adopt for the duration in Gyr of the star formation 
burst Tburst = tgso — t* the following parametrization; 

^ in*, ifM,ph>M,v 
n^MM.,,) - j ^ Mj^^ .^^^^^^ ^ ^.^^ (8) 

where = 0.5 Gyr is the burst timescale for galaxies with 
Msph > M^j,^ {M*pf^ ~ 1.5 X 10" Mq is the mass of an 
L% ~ 3. X 10^° 1/0 elliptical galaxy). In the following we 
will refer to this parameterization as case A. We also tested 
a case B, in which star formation lasts Tburst =0.5 Gyr for 
all the spheroids, independently of the mass. It is worth 
noticing that Monaco et al (2000), in order to account for 
the observed statistics of QSOs and elliptical galaxies in the 
framework of hierachical structure formation, introduced a 
a time delay decreasing with mass between the beginning of 
the star formation and the QSO bright phase. 

It is possible to explore whether the rates of star forma- 
tion implied by equation ^ are physically plausible in proto- 
galaxies. In a dark halo with circular velocity Vc and mass 
Mh oc V^, the gas is fed to the star forming regions at a 
rate which is the minimum between the infall rate and the 
cooling rate (White and Frenk 1991). The dynamical time 
is Tff{r) = (3n /32Gp{r)y^^ . The cooling time is defined as 
the ratio between the specific thermal content of the gas and 
the cooling rate per unit volume 

3p,{r) kT 



(9) 



2 nTUp ne(r)2A(r) 

where A(r) is the cooling function, depending also on chem- 
ical abundance (Sutherland and Dopita, 1993), and T is the 
gas temperature. The latter is given by kT = 1/2 prUpV^ , 
where p is the mean molecular weight of the gas and Vc is 
the circular velocity of the DM halo, whose mass is Mh- Us- 
ing the relationship between the Vc, the virialization redshift 
z-air and Mh (White and Frenk 1991) we have 

2/3 



r = 35.9 



Me 



3.3 X 105 M( 



(1 



K 



(10) 



We use n appropriate for a fully ionized gas with 25% he- 
lium in mass and adopt in the following Zvir = 3 as a refer- 
ence value. The total mass of baryons is Mtot ~ fbMn- We 
adopted for the baryon fraction fb — 0.13. 



In the following we adopt a density profile for the DM 
appropriate for a ACDM cosmology (Navarro, Frenk and 
White 1997). Depending on mass and density profile, we 
define a radius Vcooi inside which the gas can cool and col- 
lapse on timescales Tburst- We denote by fcooi the fraction 
of baryons within this radius. 

In order to explore the SFR, we use the following equa- 
tion for the evolution of the mass of cool gas Mcg{t) 



Mc^it) = -{1 - R) 



M^-M„ + CeM-—)- (11) 

The first term on the r.h.s. describes the cool gas converted 
into stars on a timescale r* and the restitution of gas from 
stars. We assume that r* = max ((tc), (t//)), where the two 
times are computed using the average density inside Vcooi, 
the limiting radius of the star forming gas. 

The SFR=Mcg/r* is expected to be affected by SN feed- 
back. However the fraction of energy transferred to the gas 
by the SNae is rather uncertain. The second term in equa- 
tion |l^ is a plausible description of this feedback (Kauff- 
mann, Guiderdoni and White 1994), and represents the rate 
of production of 'warm' gas that can not form stars imme- 
diately 



Mw(t) 



4 lyfhEsNSFR 

5 Vc^ 



(12) 



where v is number of SN produced per unit mass of formed 
stars. We adopt here and in the following a Salpeter initial 
mass function (IMF), <1> oc with x = 1.35, between 



Ml 



O.15M0 and M„ 



I2OM0 and a supernova mass 



threshold Msn > 8 Mq- With these assumptions u ~ 8.7 x 
10^'' Mq^ . We set the fraction fh of the SN energy output 
(EsN — 10^^ ergs) that is transferred to the gas to 0.15; of 
course the results for low masses are rather sensitive to the 
adopted value of fh - 

The last term in equation ^ describes the rate at 
which cold gas is supplied, namely the 'infall' rate. As for 
the timescale in the exponential law, we use Tin f{r cool) = 
vsiSLK{Tc{rcooi),Tf f{rcooi))- The constant C = Mtot/nnf is 
the initial rate. 

In Fig. 1^ the SFR, estimated using equation Fll are 
presented for spheroids with halo mass Mh ~ 1.5 x l(r^ Mq 
and Mh = 1 x lO'^^ Mq . The final masses in stars are Mgph = 
4.3 X 10" Mq and Msph = 5 x 10*' Mq respectively. It 
is apparent that the SFR can be as high as 1000 Mq in 
most massive spheroids. In these cases the fraction of the 
baryons involved in the star forming process is fcooi — 0.2, 
since the gas in the outer regions has longer cooling time. At 
small masses fcooi ~ 1, but the effect of SN feedback keeps 
the SFR low and prolongues the star formation, despite the 
shorter cooling times. 

We exploited equation 11 to compute the SFR for ob- 
jects with final mass in stars ranging from 10^ to 5 x 10^^ 
Mq - The SFR, averaged over the time Tburst , as function of 
mass turns out to be 



SFR'. 



100 (AIsph/1 X 1O"M0)' ' Mq yr-\ 



(13) 



Only about 20-40% of the baryons are locked in stars, 
the rest being in form of warm gas. For large hosts only a 
small fraction of the gas has a cooling time short enough 
to collapse before the QSO reionizes, reheats and possibly 
expels the gas, stopping both star formation and accretion 
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Figure 2. Star Formation Rate for spheroidal galaxies with halo 
mass M[{ = 1.5 X 10"''' Mq and final mass in stars M^ph = 
4.3 X 10^1 Mq (upper curve) and with Mh = 1 x lO^^ Mq and 
^sph = 5 x 10^ Mq (lower curve). The change from solid to 
dashed line marks the end of the star burst at Tj^^st given by 
equation ^, because of the QSO shine. In the paper we assume 
that the SFR vanishes for T > Ti,urat ■ 



on the black hole. For small mass hosts, star formation can- 
not proceed much before the QSO shining because of the 
efficiency of SN feedback. As a result the stars are again a 
small fraction about 3-5% of the DM mass. 

We checked that the average stellar metallicity increases 
with the mass. Large SFR in massive hosts rapidly increases 
the metallicity of the gas and forming stars in such a way 
that it is about solar before the QSO advent. Conversely 
in small hosts, the SNe feedback keeps the average SFR 
low before the QSO blocks the star formation. Thus, al- 
though the QSO advent is delayed, nevertheless a lower 
metal fraction is produced. In our scenario SNIa explode 
in large hosts after the star formation and QSO output pre- 
vent the gas from forming stars. For less massive objects 
still important episodes of star formation can occur even 
after the SNIa begin to explode. This produces a positive 
E/Fe-(T correlation, as observed. The slope of the correlation 
significantly depends on SNIa progenitors and the ensuing 
timescale for Fe enrichement (Thomas, Greggio and Bender 
1999; Kobayashi, Tsujimoto and Nomoto 1999; Nomoto et 
al 1999). The details of the chemical evolution will be dis- 
cussed elsewhere (Romano et al. in preparation). 

The rate at which spheroids of present day mass Msph 
begin to form stars at is directly related to the rate at 
which they appear to host a QSO (equation 6 and Fig. |l]) 
at time tQso = + Tburat{M). 



3 THE SPECTRAL EVOLUTION OF 

SPHEROIDAL GALAXIES AND THEIR 
FAR-IR EMISSION 

In the previous section we computed the cosmological rate 
at which spheroids begin their star formation activity. We 
also estimated how the SFR depend on M^ph- To compare 
with observations we need as well predictions about their 
SED evolution. 

The galaxy surveys at submillimeter wavelengths with 
SCUBA and analogies with local starburst galaxies suggest 
that the starburst phase occurred in a dusty environment. 
The spectrophotometric evolution in the UV and optical 
bands of bursts of star formation, under the assumption of 
no or negligible dust absorption, has been investigated by 
many authors (e.g. Bruzual 1983; Arimoto & Yoshii 1986; 
Bruzual & Chariot 1993; Bressan, Chiosi & Fagotto 1994). 
While the no dust approximation well describes the physical 
situation in elliptical galaxies at later times, dust is expected 
to play an important role during the starburst epoch. The 
main effect of dust is to transfer a major fraction of the 
emitted power from the UV and optical bands to the mid- 
and far-IR. 

In order to compute the SEDs of the spheroids with vig- 
orous star formation and dust, we use the model GRASIL 
developed by Silva et al. (1998), which is state-of-the-art 
as far as dust reprocessing is concerned. GRASIL has been 
tested against UV to radio SEDs of local spirals and star- 
burst galaxies (Silva et al. 1998, Silva 1999). It has also 
been used to reproduce the IR properties of local galaxies 
in the framework of semi-analytical models (Granato et al 
2000). GRASIL includes: (i) chemical evolution; (li) dust 
formation, assumed to follow the chemistry of the gas; [Hi) 
integrated spectra of simple stellar populations (SSP) with 
the appropriate chemical composition; and [iv) appropriate 
distribution of stars, molecular clouds (in which stars form 
and subsequently escape) and difi^use dust. In particular, 
dust is spread over the whole galaxy and its temperature 
distribution is determined by the local radiation field. It is 
noteworthy that the assumption of a single dust tempera- 
ture turns out to be a significant oversimplification. 

The predictions of sub-mm fluxes are plagued by the un- 
certainty of the far-IR wavelength dependence of the dust 
emissivity which, convolved with the temperature distribu- 
tion of dust grains, determines the sub-mm steep slope of 
the SED. Longward of ~ 40/im the dust opacity can be rep- 
resented by a power law fc^ oc Values ranging between 
1 and 2 are commonly adopted for the spectral index 7. 
The classical computations by Draine and Lee (1984) pre- 
dict 7 = 2, but there are several indications of a shallower 
slope at least in some cases. The issue is further confused by 
the fact that in most papers the observed sub-mm contin- 
uum is fitted with single temperature gray bodies. However, 
Silva et al. (1998) found that the IR spectrum of the archety- 
pal ultraluminous infrared galaxy (ULIRG) Arp 220 is well 
represented by k^, cc v^'^ , whilst oc is more suitable 
for M82, the prototype starburst galaxy, and other galaxies. 
Arp 220 may however be more representative of a spheroid. 
Since a single opacity law may be too strong an assumption 
we treat 7 as a free parameter allowed to vary between 1 
and 2. 

In Figs. ^ and ^ model SEDs for elliptical galaxies 
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Figure 3. Model SEDs of elliptical galaxies before (thicker lines) 
and after (thinner lines) the QSO activity for Ti^urat = 0.5 Gyr. 
Numbers along the curves are the age in Gyr of the corresponding 
models. 




log X[fim] 

Figure 4. Same as the previous figure but for T),y^rst = 2 Gyr 

are shown as a function of time adopting 7 = 1.5. In the 
former case, gas and dust are removed and star formation 
is stopped at Tburst ~ 0.5 Gyr (correspinding to the SFR in 
the upper curve of Fig. ^ , while in Fig. ^ Tburst = 2 Gyr 
(lower curve of Fig. ^). It is worth noticing that GRASIL 
predicts that the fraction of bolometric luminosity emitted 
in the far-IR decreases with time during the burst phase. In 
the case depicted in Fig. ^this fraction declines from about 
90% to about 40%. 

At the end of the burst, gas and dust are removed and 
the stars already formed emerge and shine in the optical 
bands, while the far-IR luminosity drops dramatically. Of 
course both cases exhibit the same behavior at r Tburst- 
However, the model with Tburst ~ 0.5 Gyr passes through 




Figure 5. Predicted contributions of the spheroids to 850 ^m 
counts. Solid line is for case A {Tburst given by equation H) and 
ki^ cx u^-^) while dot-dashed line is for case B {Tburst = 0-5 Gyr 
and fci/ oc u^'^). Data are from Blain et al. (1999b) (open circles), 
Hughes et al. (1998) (star and open triangles), Barger, Cowie and 
Sanders (1999) (open squares) Bales et al. (2000) (filled circles), 
Chapman et al. (2000) (filled triangle) Borys et al. (2000) (filled 
squares). 



a relatively blue phase (in the rest frame) during the first 
Gyr after the onset of the galactic wind, while in the case 
Tburst ~ 2 Gyr, the galaxy emerges with redder colours. 

The UV and optical luminosities of the spheroids during 
the starburst dusty phase critically depend on the relative 
distribution of stars and dust, which sets the small fraction 
of UV luminosity escaping from molecular clouds (where 
stars are born). Therefore any model has large uncertain- 
ties, and only observations can enlighten the physical pro- 
cesses involved. For instance, local starburst galaxies, which 
are both strong IR emitters and relatively blue objects, are 
well reproduced by GRASIL, allowing stars to emerge from 
rather thick clouds (r(l/j,m) ~ 30, as typical in the Galaxy) 
on timescales of the order of 10^ yr (Silva et al. 1998) . Con- 
versely Silva (1999) has shown that a very red optical spec- 
trum, resembling that of the Extremely Red Object HR 10 
(Hu & Ridgway 1994; Cimatti et al. 1998; Dey et al. 1999), 
can be obtained assuming longer timescales for the escape 
of stars from molecular clouds, endowed with much lower 
optical thickness (r(l/im) ~ 2). On the other hand the SED 
predictions are more robust in the far-IR, the most impor- 
tant uncertainty beeing related to the dust opacity spectral 
index 7 discussed above. 

The radio emission is computed starting from the type 
II SN rate (Condon & Yin 1990), under the assumption that 
all stars with mass Msn > 5 Mq end in a SN; the radio 
luminosity has to be decreased by a factor ~ 2 if Msn > 
8 Mq is assumed (see Silva 1999 for details). 
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Redshift 

Figure 6. Predicted redshift distribution of spheroids brighter 
than 1 mjy and 10 mjy at 850 fim for case A. 



4 THE SUBMILLIMETER COUNTS AND THE 
DURATION OF THE BURST IN SPHEROIDS 

In our scheme we identify the observed galaxies with 
"Ssso ^ 1 nijy with the dusty starbursts, occurring in 
spheroidal galaxies before the shining of the QSO. From 
Figs. ^ and^ it is apparent that the duration Tburst of the 
dusty starburst can be tested by far-IR counts and related 
statistics. 

Coupling the SED predicted by GRASIL with the for- 
mation rate of the spheroids n-i,{Msph,t*), we compute the 
expected contribution of spheroids to the submillimeter 
galaxy counts. The parameters mainly affecting the fit to 
observed counts are the spectral index of dust emissivity 7 
and Tburst- 

With the assumption Tburst ~ 0.5 for all masses, the 
spheroids would reproduce 850 ^m counts only at bright- 
est fluxes, allowing for a substantial contribution from an 
additional population of possibly lower redshift objects. By 
converse with a burst duration depending on the mass as in 
case A (equation^ we obtain an excellent fit to the observed 
counts with 7 = 1.5, as shown in Fig. |^. Therefore in the 
following we concentrate on this solution, unless otherwise 
explicitly stated. 

It is worth noticing that also a longer duration of the 
burst Tburst = 2 Gyr, with a corresponding decrease in 
SFRs, could give an acceptable fit to the counts, provided 
that 7 = 2 is adopted. Howevever this solution is not sat- 
isfactory, since it would decrease the the Mg/Fe ratio be- 
low the average observed value for local elliptical galaxies 
Mg/Fe ~ 0.2 (Kobayashi and Arimoto 1999; Trager et al. 
2000 a,b). 

Since in our scheme the infrared emission in the 
spheroids occurs before the QSO shine, the SCUBA sources 
must be located at high redshift. For fiuxes Sgso ~ 1 mJy, 
the predicted redshift distribution is relatively broad, with 
almost all the sources lying in the redshift interval 2 ^ z ^ 5 

(Fig. However we predict a non negligible number of 



sources at z ^ 5. Even at a fiux limit Sgso ~ 10 mJy almost 
all the sources are predicted to lie at z ^ 2. This is in agree- 
ment with the results obtained so far by the spectroscopy 
of the optically identified counterparts of SCUBA sources 
(Blain et al. 1999c; Small et al. 2000). A testable prediction 
of the proposed scheme is that there are more extragalactic 
sources with 5*850 Si 1 mJy at a > 5 than at z < 1, where 

only 1% are located. 

The model source counts at Sgso ~ 1 uijy are domi- 
nated by the massive spheroids (E/SO galaxies), which have 
presently large luminosities Mb < Mb* — —20.5. 

How do the properties of the starbursting spheroids 
compare with those of the individual submillimeter galax- 
ies? The estimates of SFR for SCUBA galaxies range from 
several hundreds to a few thousands M0yr~^ (Ivinson et al 
2000; Frayer et al 2000) Two SCUBA sources have been op- 
tically identified with Extremely Red Objects. The available 
optical identifications point toward very faint (7 > 26 — 27) 
counterparts (Smail et al. 1999). In our model, the most 
massive spheroids have SFR ranging from several hundred 
to thousand Mqjt'^ and they are faint at UV and op- 
tical wavelengths during the dusty phase. However, it is 
worth mentioning that the expected optical magnitudes and 
colours of these high z objects are very sensitive to the 
tiny fraction of UV emission surviving after dust extinc- 
tion and can in practice span a large range of values, even 
at a fixed redshift. For instance, we can reproduce the 850 
^m source counts with objects exhibiting, during the dusty 
phase, colours similar to those of EROs, as well as with rela- 
tively bluer ones. This can be achieved simply by tuning the 
small fraction of young stars which have already moved out 
from the parent clouds and/or on the typical optical depth 
of clouds (cfr. sect. 3 and 5). 

The extraordinary bursts of star formation before Iqso 
result also in large radio luminosities. Our scenario predicts 
that the 850 pm sources with 1 ^ 5*850 ~ 10 mJy have radio 
fluxes at 5 GHz in the range 5 ^ S5 ^ 200 /xJy, in agreement 
with the first results on the radio emission of SCUBA galax- 
ies obtained by Smail et al. (2000) and with the findings of 
Barger, Cowie and Richards (2000). 

At 450 Blain et al. (1999c) report N{> S = 
10 mJy) = 2.1±1.2x 10^ deg"^. We predict that at this flux 
level the dusty spheroids are 0.7 x 10'' deg~^, a significant 
fraction of the counts. By contrast, at 175 ^m the expected 
contribution is negligible at the bright limit 5i75 ~ 100 mJy 
(Puget et al. 1999; Kawara et al. 1998). This is due to the 
fact that the K-correction is not as favourable as it is at 
longer wavelengths, dimming sources at z Si 3, where most 

of the dusty spheroids are located (cfr. Fig. Indeed Scott 
et al. (2000) have observed with SCUBA at 450 and 850 fim 
a subsample of the objects selected at 175 ^m. The result- 
ing submillimeter spectral energy distributions suggests that 
these objects are at low or moderate redshift < z < 1.5. 

Going to shorter wavelengths, the spheroids emerge as 
an important component below 10 mJy at 60 /xm; at 90 /im, 
they start to be a non negligible fraction at fluxes ^ 20 mJy. 
The latter limit is within the reach of the ELAIS survey, 
covering about 20 square degrees at 90 fim (Efstathiou et 
al. 2000). 

Several authors presented models of galaxy evolution 
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taylored to fit available data and in particular the submil- 
limeter counts. For instance Guiderdoni et al (1998) repro- 
duced the counts by imposing a large evolution in the frac- 
tion of mass density involved in bursts of star formation 
and in the fraction of Ultra Luminous Infrared Galaxies 
(ULIGs), both fractions passing from a few percent in the 
local Universe to substantial values at high redshift. The 
conclusion that rapid star formation activity is required at 
z ^ 3 has been claimed by Blain et al. (1999a). They showed 
that, in order to reproduce the observed sub-mm counts, 
the activity in galaxies had to be a factor of 200 larger at 
z ^ 3 than in the local Universe. In their scheme this activ- 
ity refers to both starbursts and to AGNs, and is triggered 
by merger events at high redshift in a hierarchical model of 
galaxy formation. A view even closer to our own has been 
proposed by Tan, Silk & Balland (1999). They exploited a 
model in which the final morphology of the galaxies depends 
on the number of collisions and tidal interactions they suf- 
fer (Balland, Silk & Schaeffer 1998), deriving for spheroids 
large star formation rate density at high redshift. In particu- 
lar, they broadly reproduce the 850 fim counts with objects 
that suffered important collisions and that they identify as 
precursors of spheroidal objects. 

Our model is characterized by a very steep increase of 
the counts below ~ 60 mjy. Above this flux limit the sources 
are mainly low redshift spirals and starburst galaxies, while 
just below it a large amount of high redshift star forming 
spheroids appear. 

In conclusion the most massive spheroids bursting at 
2^2, before the onset of the QSO activity, can be iden- 
tifled with the observed SCUBA sources down to 1 mJy. 
The SCUBA galaxies are strongly clustered, since they are 
shining in the most massive virialized haloes at high red- 
shift. The large rate at which the spheroids form at high 
redshift, is dictated by the QSO luminosity function and by 
the short QSO lifetime AtQ . The other relevant parameter is 
the duration of the star formation, observationally bounded 
also by the chemical evolution of QSO hosts and in general 
of elliptical galaxies. The implied star formation rates are 
physically plausible in protogalaxies (cfr Sect. 2.2). 



5 STARBURSTING AND POST-STARBURST 
SPHEROIDS AND LYMAN BREAK 
GALAXIES 

In the context of the star formation in spheroids at high red- 
shift, a relevant issue is the relation of high redshift galaxies 
selected in optical and near IR bands with the submillime- 
ter selected galaxies. Steidel et al. (1999), using a galaxy 
sample complete to Iab ^ 25 and appropriate color crite- 
ria, were able to identify a large sample of galaxies with the 
Lyman Break falling in the optical or near UV bands. These 
objects at 2; > 3 have been claimed to exhibit dust absorp- 
tion with reddening in the range < E{B ~ V) < 0.4, with 
an average value of ~ 0.15 (Steidel et al. 1999). After cor- 
rection for the implied dust extinction, the UV luminosity 
density of the LBGs at 2 ~ 3 — 4 would be almost equal to 
that inferred from galaxies lying at 2 ~ 1 (Connolly et al. 
1997). The corresponding extinction corrected (by a factor 
of about 5) star formation rate of a typical galaxy would 




Figure 7. Track of the model of Fig. 1, assumed to begin star 
formation at 2* = 5, in the two colour plane (G — 1Z) vs. {TZ — I) 
used by Steidel et al. (1999) to select LBG candidates. The upper 
left region bounded by solid lines is the selection region. Number 
along the track (solid line before the wind and dotted afterwards) 
are the corresponding redshifts. 



be of about 60 Meyr"^ . This implies a SFR per unit volume 
almost constant psFR — 0.15 Mq yr~^ Mpc~^ from 2 ~ 1 
to 2 ~ 4. 

On the other hand SCUBA observations of a selected 
subsample of LBGs resulted in only one possible detection, 
suggesting that the SFR of the population as a whole is 
possibly low (Chapman et al. 1999). However this result is 
expected in our scheme, because of selection effects (see also 
Adelberger and Steidel 2000; Peacock et al. 2000). The UV 
luminosity during a star burst critically depends on the frac- 
tion of young stars already escaped from their parent molec- 
ular clouds and on the associate time scale (cfr. Section 3 
and Silva et al. 1998). We already remarked that it is more 
prone to model uncertaities than the IR and sub-mm emis- 
sion. In any case, the typical behavior of our models is shown 
in Fig. |l] and Fig. The ratio Lfir/ Luv varies from 5-10 
at about the end of the burst to very large values just after 
the begin. After the QSO phase the ratio falls down very 
rapidly, in the short period when the spheroids are still blue 
unextincted objects. 

Keeping in mind the above mentioned caveats concern- 
ing the dusty star-forming phase, our modelling predicts 
that the spheroids may spend most of their 2 > 3 life inside 
the region of the colours plane used to select the LBG can- 
didates (Steidel et al 1999). For instance, a massive galaxy 
Msph — 3 X 10^^ Mq with star formation started at 2* ~ 5 
and with the QSO activity appearing after Thurst — 0.5 Gyr, 
would have Iab ^ 25 since 2^4, and would exhibit colours 

(G — 7V) and {TZ — X) quite similar to those of the LBGs, as 
evidenced by Fig. ^ 

In our scenario the hosts of bright QSOs (Lhoi > 
10^"^ 1/0 may appear as LBGs both during their starburst 
phase and soon after the QSO shining; these hosts now have 
luminosities Mb — 20.5 and velocity dispersion a ~ 200 
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km/s. As stated in section 4, these hosts are also seen as 
SCUBA objects during their starburst, since they have av- 
erage a SFR in the range 100-1000 M© yr"\ The metal 
abundance of their stars, estimated with a Salpeter IMF, 
is about solar. These massive spheroids are also expected to 
exhibit a large rate of SNIa explosion, after the QSO shining, 
when they are almost gas and dust free. 

Less massive spheroids, with present luminosities — 17 < 
Mb < —20.5, at redshift z ~ 3 have average SFR from a few 
to one hundred Mq yr^^ and sub-solar metal abundance 
in stars, decreasing with decreasing SFR. The majority of 
these spheroids have not yet experienced the QSO shining at 
z > 3 and they are in the pre-QSO phase. In our modelling 
of the SED (cfr. sect 3) dust follows the gas chemistry and, 
as a consequence, the extinction is expected to decrease with 
decreasing SFR. Thus the less massive spheroids can appear 
as LBGs and may constitute a major fraction of them. In 
conclusion the spheroids independently of their mass very 
likely showed up at 2; > 3 as LBGs. However they have 
SFR, metal abundance, dust extinction and clustering scale 
decreasing with their mass. 

It is also worth noticing that in our schematic model 
star formation has completely stopped after tqso- However, 
ionization, heating and winds due to the QSO can simply 
decrease the SFR from several hundred Mq yr~^ to several 
Mq yr~^, compatible with the observed UV emission from 
a typical LBG. Since this more prolonged star formation 
would involve only a minor fraction of the galaxy mass, say 
5-10 per cent of the mass in 1-2 Gyr, the present day spectra 
of the galaxy would be unaffected. A residual low level of 
star formation activity has been inferred from broad band 
spectra of the HDF elliptical gala:xies at « ~ 1 (Franceschini 
et al. 1998). 



6 THE QSO PHASE 

The QSO can heat, ionize and also expel the ISM, thus 
strongly decreasing or even stopping the copious star for- 
mation in the host galaxy at tqso (Silk and Rees 1998; 
Fabian 1999). However still significant star formation can be 
present during the QSO lifetime Afg ~ 4 x 10^ yr. Several 
objects with these characteristics have already been found 
(e.g. Omont et al. 1996; Ivison et al. 1998; Benford et al. 
1999; Yun et al. 1999). If we assume that the star formation 
is not quenched during the nuclear activity, then the fraction 
of QSOs in an unbiased submillimeter survey depends on 
the ratio between the QSO lifetime and the burst duration 
Tburst- Since the latter timescale in our scheme is a function 
of the mass, this precise fraction depends on the flux limit. In 
particular for bright submillimeter sources Sgso ~ 2 mjy we 
expect Tburst — 0.5 — 1 Gyr and then a fraction of QSOs from 
a few to 10% is estimated. Recent observations of SCUBA 
galaxies with Chandra suggest that they are mainly powered 
by starbursts, unless Compton-thick tori with little circum- 
nuclear X-ray scattering are in place (Fabian et al 2000; 
Hornschemaier et al 2000; Severgnini et al 2000). A rather 
different result has been obtained by Bautz et al (2000), 
who were able to detect the X-ray counterparts of two sub- 
millimeter sources in A 370. These authors conclude that 
20^jg% of the submillimeter galaxies exhibit X-ray emis- 
sion from AGN. Spectroscopy revealed no QSO among the 



15 850 fim selected sources, though several of them exhibit 
possibile signs of low nuclear activity (Ivison et al. 2000). 

On the other hand hard X-ray selected sources down to 
S2-iokeV ~ 2 — 4x 10~^^ ergs s~^ cm~^ are basically not de- 
tected even in deep submillimeter exposures Sgso ~ 2 mJy 
(Barger et al 2000; Fabian et al 2000; Severgnini et al 2000), 
with a few exceptions. In particular in the sample of Barger 
et al (2000) most of the X-ray contribution to the resolved 
X-ray background comes from AGN with Uoi ~ 10""^ - 10'"' 
erg s~ , associated to BH masses Mbh ~ 10^ - 10* Mq and 
baryon masses in spheroids Mgph ~ 5x10^— 5x10^" Mq , in 
keeping with the predictions of Salucci et al (1999). At 850 
Hm, contributions from dusty tori around the active nucleus 
and from star bursts in host galaxies are expected. However 
the dust emission from the torus peaks at A ~ 30 ^m, though 
in case of very extended structures a wide spread of temper- 
atures tends to produce a plateau up to 80 fim (Granato and 
Danese 1994; Efstathiou & Rowan- Robinson 1995). Thus 
relatively faint AGN at z ~ 1 — 2 can not be detected at 
Ss50 > 1 mJy. Also the starburst emission is elusive, be- 
cause of the relatively low star formation expected in these 
objects. For instance in our model we expect that their host 
galaxies have already formed most of their stars in about 
1-2 Gyr, with an average SFR ~ 2.5 - 50 Mq yr~^ . These 
limits are below possible detection by SCUBA for galaxies 
at z > 1. Then, as claimed by Granato, Danese & Frances- 
chini (1997), the obscured AGNs, which produce a major 
fraction of the Hard X-ray Background Radiation (HXRB), 
contribute only a small fraction of the FIRB at submillime- 
ter wavelengths (see also Gunn ands Shanks 1999). 

Interestingly enough, the evolutionary links between the 
star bursting phase and the nuclear activity we have envis- 
aged at high redshift, can be applied to ultra luminous in- 
frared galaxies at low redshift (Sanders et al. 1988). In our 
scheme the low redshift ULIRGs correspond to the dusty 
starburst phase. This is in keeping with recent ISO observa- 
tions that show that for the local ULIRGs, the starburst ac- 
tivity is sufficient to power the bulk of the bolometric lumi- 
nosity (Genzel et al. 1998). However we also expect cases in 
which the QSO light reprocessed by dust in the host galaxy 
is dominating over the far-IR emission due to the starburst. 



7 DISCUSSION 
7.1 QSO hosts 

A major assumption in our model is that the relationship be- 
tween BH mass and host galaxy mass has been imprinted at 
early epochs (section 2.1). Recent observations from ground 
and with HST show evidence that the radio loud QSO and 
the radio galaxies at high redshift z > 1 are hosted in bright 
spheroidal galaxies with a rather evolved stellar population 
(Dunlop 1999; Zirm et al 1999; Lacy, Bunker and Ridgway 
2000; McLure and Dunlop 2000). These observations sup- 
port our assumption. Also recent photometric and spectro- 
scopic follow-up of a hard X-ray selected sample showed that 
a large fraction of the z ^ 1 AGN, responsible for a major 

fraction of the HXRB, are hosted in early type optically 
luminous galaxies (Barger et al 2000). 

Observations of radio quiet QSO hosts show that their 
luminosity at fixed QSO luminosity decrease with increasing 
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z, (McLure at al 1999; Ridgway et al 1999; Rix et al 1999). 
The 5 QSOs with z > 1.7 studied by Rix et al (1999) exhibit 
a median ratio of the nuclear to the host galaxy luminosity 
LN/Lhost — 6 in the V-band, while in the R-band for 9 
low redshift {z < 0.24) QSOs McLure et al (1999) found 
LM/Lhost — 1.5. This decrease by a factor 4 is explained 
by the observed decline of the median Eddington ratio from 
L/LEdd ~ 1 - 2 to L/LEdd ^ 0.1 - 0.2 for QSO at high red- 
shift and low redshift respectively (Padovani 1989; Sun and 
Malkan 1989; Wandel 1999). As a matter of fact the low red- 
shift radio quiet QSO of the sample of McLure et al (1999) 
exhibit a median L/LEdd — 0.2, under the assumption that 
Mbh/Msph 3 X 10"^ 

The observed decrease of the L/Lsdd ratio is expected 
in our scenario. When QSOs shine at high redshift z > 1.5 — 2 
they have still large gas reservoirs in the hosts and they 
can reach the self-regulated Eddington accretion. At lower 
redshift, gas in hosts is scanty and only encounters with 
other galaxies can supply material for re-activation. However 
this supply can only occasionally be large enough to reach 
the self-regulated regime (see also Cavaliere and Vittorini 
2000). Note that the decrease of L/Lsdd is quantitatively 
reproduced by the semi-analytic model by Kauffmann and 
Haehnelt (2000). In this case however it is driven not only 
by the lower gas fraction left at lower redshift, but also by 
the adopted redshift dependence of the accretion timescale 
tacc oc (1 -I- z)~^'^ in their best fit model. 



7.2 Early baryon collapse in massive DM haloes 

The hierarchical scheme in the Dark Matter component can 
naturally explain the onset of QSO activity and of star for- 
mation in spheroids. In particular the Press-Schechter for- 
malism within the standard variants predicts a large enough 
number of massive DM halos in order to host QSOs and 
their parent galaxies (see e.g. Haehnaelt and Rees 1993; 
Haehnaelt, Natarajan and Rees 1998; Cattaneo, Haehnaelt 
and Rees 1999; Monaco, Salucci and Danese 2000). However 
a distinct character of hierachical DM scheme is that smaller 
haloes virialize earlier. On the other hand the high SFRs and 
IR luminosities implied by the far-IR counts exclude that 
most of the stars have formed in smaller subclumps at very 
high redshifts, say z ~ 10. It is plausible that at high redshift 
z ^ 6 the SFR has been kept very low in the subclumps by 

'internal' mechanisms (e.g. SN feedback or even a low activ- 
ity of the nucleus), which became progressively ineffective 
when larger and larger masses virialized. Then, as discussed 
in Sect. 2.2, a significant fraction of the baryons in massive 
'virialized' haloes can fragment very rapidly to form stars 
and QSOs, while the star formation is further delayed in 
smaller haloes by the effect of stellar feedback. The shining 
of the QSOs is the final and brief time episode of the complex 
processes involved in assembling baryons in stars and BH. 
Most of the relevant observational facts can be framed in a 
scenario in which the most massive objects complete their 
activity in a very short time interval Ti,arst — 0.5 Gyr and 
the less massive in ~ 2 Gyr. In a sense the DM hierarchi- 
cal sequence is corrected by an Anti-hierarchical Baryonic 
Collapse (see also paper II), regulated by the QSO and SN 
feedback. 

Since star formation in massive ellipticals occurs early 
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Figure 8. Star formation (solid line and data points, left scale) 
and mass accretion (dotted line, right scale) rates per unit volume 
as functions of redshift. The latter has been derived from the data 
on QSO luminosity function evolution as a function redshift. The 
SFR refers to the Case A. 



and in relatively short times, the emerging picture is some- 
what in between the concept of a 'monolithic' collapse of 
baryons inside the haloes and the 'classical' hierarchical 
semi-analytic models, discussed for instance by Kauffmann 
& Chariot (1998) and Governato et al. (1998). In particular, 
the connection between ellipticals and quasars has been ad- 
dressed recently by Kauffmann & Haehnelt (2000). In their 
'classical' scenario of the merging process they reproduce the 
QSO activity decline, but the most massive spheroids form 
late. As a result they predict a significant active evolution 
of the host galaxy after the quasar shines. On the contrary 
in our scheme the host galaxies evolve passively, except for 
subsequent minor mergers at lower redshift which change 
the mass in stars by no more than 20% (Trager et al. 2000b; 
see also Van Dokkum et al. 1999) . 

It is very likely that the main aspects of the Anti- 
hierarchical Baryonic Collapse can be reproduced by semi- 
analitycal models, by acting on very basic parameters, such 
as merging and star formation timescales and SN feedback 
efficiency. 



7.3 The star formation in spheroids 

The star formation rate per unit volume in spheroids as a 
function of redshift derived from our modelling is presented 
in Fig. ^. The mass accretion rate per unit volume implied by 
the observed nuclear activity is also presented. It is apparent 
that in our scheme the star formation rate in QSO hosts is 
rapidly falling down at z ^ 2. However the star formation 

rate in the Universe remains significant because of the star 
formation in discs and dwarves. 

The star formation rate inferred, without dust correc- 
tion from the UV emission of LBGs galaxies at z > 3 is a 
factor ~ 5 lower than our prediction (Steidel et al 1999). 
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However this is consistent with the evidence of strong ex- 
tinction of UV emission. 

The density of mass cycled through stars in QSO hosts 
amounts to Pc* — 3.2 x 10* Mq Mpc~^, which implies 
(using Salpeter IMF) a present density of long-lived stars 
in QSO hosts p* ~ 2 x 10* Mq Mpc"^. This mass den- 
sity can be compared to the mass density in stars in local 
spheroids 2 x 10* < < 5 x 10* Mq Mpc"^ (Fukugita 
et al 1999). The metal mass density is roughly given by 
pz — yzpc* — 9 x lO'' Mq Mpc~'^, yz denoting the fraction 
of metals returned to the interstellar matter by a stellar 
generation (0.028 with the adopted IMF). In our scenario 
only about 40% of these metals are still locked into stars, 
while the rest are spread over the intergalactic medium. 
The X-ray emitting gas in groups and clusters of galaxies 
amounts to Pigm{,z = 0) ~ 1 x Mq Mpc~^ (Fukugita 
et al. 1999). Available observations suggest that this gas 
has Z ~ I/ZZq (Renzini 1999), implying pz igm{z = 0) ~ 
7 X 10*^ Mq Mpc"^, not far from our prediction. 



8 SUMMARY 

We have presented a unifying scheme for the formation of 
QSOs, elliptical and SO galaxies, and for the bulges of spirals. 
We have shown that the main aspects of the evolution of 
spheroidal galaxies and QSOs can be well understood if we 
assume at t = tqso, when the QSOs shine, the bulk of 
the stars have already been formed in QSO hosts with star 
formation lasting a time Tburst ~ tqso — t* ranging from 
~ 0.5 to ~ 2 Gyr when going from more to less massive 
hosts. This suggest an Anti-hierarchical Baryonic Collapse, 
i.e. that the baryons in large spheroids collapse very rapidly 
to form stars and QSOs, while the collapse of baryons in 
smaller spheroids is slowed down. In fact we have shown 
that the central regions of the massive DM haloes can form 
stars on very short timescales, while the SNe feedback is 
strongly slowing down the SFR in small spheroids, extending 
the star forming phase. The ensuing formation rates scale as 
SFR ~ 100 (M/IO^M©)^'^ Mq yr'V In our picture when 
the QSO reaches its maximum, the combined actions of its 
power and SN feedback are able to practically stop both the 
SFR and the nuclear accretion through ionization, heating 
and winds. This influences both the mass in stars Msph and 
the BH mass Mbh and may be the origin of the observed 
correlation between the two masses. As a result only about 
20-40% of the baryons are presently in stars both in massive 
and small DM galactic haloes. The remaining is in warm or 
hot gas with Z ~ \Zq. The outlined scenario can explain: 

(1) the main aspects of the chemical evolution of 
spheroids (stellar metallicity, luminosity-metallicity rela- 
tionship, the a enhancement) and the observed elemental 
abundances in QSOs; 

(2) the 850 source counts; 

(3) the main properties of Lyman Break Galaxies; 

(4) the detection of high redshift, red and old massive 
ellipticals. 

This general view implies the following evolution se- 
quence for QSO hosts: 

(i) the most massive appear, at high redshifts z ^ 2 — 6, as 



ultraluminous far-IR galaxies (SCUBA galaxies). This phase 
lasts for Tburst ~ 0.5 — 1 Gyr. On the average they are 
optically faint Iab 26 — 27 with colors difiicult to predict, 

in the sense that, depending on uncertain details of dust 
distribution, both objects with some intrinsic UV emission 
as LBGs and very red objects (like EROs) can be expected; 

(ii) the hosts of intermediate and low mass have bursts 
more protracted Tturst ~ 1 — 2 Gyr and at z > 3 can appear 
as LBGs; 

(iii) the QSO phase then follows. The IR emission is at- 
tributable both to the starburst and the nuclear activity. 
This phase lasts for the QSO duty cycle Atq ~ tduty — 
4 X lO'^ yr. It ends because heating, and ionization of the 
ISM and the possible onset of galactic winds drastically re- 
duce SFR and nuclear activity; 

(iv) a long epoch of passive evolution follows, with spectra 
becoming quickly red. 

The scheme implies that 850 pm selected galaxies with 
5850 > 1 mjy are strongly clustered. Also LBGs are ex- 
pected to be significantly clustered, though to a lower level. 

The emerging picture is consistent with hierarchical 
models for structure formation. However it evidences that, 
in order to understand key observations of QSOs and ellip- 
tical galaxies, the hierarchical assembly of DM haloes must 
be completed by the history of the baryons inside the haloes 
themselves. The physical processes related to star formation 
and feedback suggest that baryons can collapse in stars and 
QSO more rapidly in more massive haloes. However why the 
QSOs shine in their hosts when star formation has already 
progressed, is still an open problem. 
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